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ABSTRACT 

The effects of seeing on the parameters of the Sersic profile are studied in an analytical 
form using a Gaussian point spread function. The surface brightness of Sersic profiles 
is proportional (in magnitudes) to r^'". The parameter n serves to classify the type of 
profile and is related to the central luminosity concentration. It is the parameter most 
affected by seeing; furthermore, the value of n that can be measured is always smaller 
than the real one. It is shown that the luminosity density of the Sersic profile with n 
less than 0.5 has a central depression, which is physically unlikely. Also, the intrinsic 
ellipticity of the sources has been taken into account and we show that the parameters 
are dependent when the effects of seeing are non- negligible. Finally, a prescription for 
correcting raw effective radii, central intensities and n parameters is given. 

Key words: atmospheric effects ~ methods: data analysis - galaxies: distances and 
redshifts - galaxies: photometry. 



1 INTRODUCTION 



Galaxies exhibit many different morphologies and therefore 
dynamical properties. They also vary enormously in size 
and luminosity. Galaxies can consist of one or more dynam- 
ical substructures. The classification methods are mainly 
ll" ' based on the morphologies of these substructures. The most 
. 5^ I widely used classification scheme is that introduced by Hub- 
^^ • ble (1936), which is based on the ratio of the spheroidal bulge 
5_j ' and disc luminosities. The galactic morphological types vary 
Cu ^ from ellipticals (with only spheroidal components) to late- 
type spirals, with small spheroids and prominent disc com- 
ponents. Irregular galaxies, which are later than spirals in 
the Hubble sequence, are characterized by the absence of 
symmetry. In recent decades much work has been done to 
improve the method for determining the principal compo- 
nents of galaxies (e.g. Prieto et al. 2000 and references 
therein). 

It is usually assumed that the light of a galaxy follows 
the mass distribution. The mass distribution can then be 
inferred by modelling the light distribution. The light of a 
galaxy is usually modelled by fitting the surface brightness 
profile of each structural component with certain analyt- 
ical curves. These laws have certain free parameters that 
must be determined during the fitting process. From de Vau- 
couleurs (1948), it is well known that the surface brightness 
profile (in magnitudes) of elliptical galaxies is proportional 
to r^'^ (where r is the radial distance to the centre of the 
galaxy) . This law was also applied to bulges of spiral galaxies 
which have similar shapes, colours and kinematics to those 
of ellipticals. It has recently been discovered that not all 



bulges follow an r^' * profile (Caon, Capaccioli & D'Onofrio 
1993; Andredakis, Peletier & Balcells 1995). A better an- 
alytical form is the Sersic profile (Sersic 1968) which has 
surface brightness proportional to r'^' " and generalizes the 
r^' ** law. The parameter n is one of the three free parame- 
ters and it defines the type of the profile. When n = 1 the 
surface brightness profile is exponential. Increasing values of 
n give more centrally concentrated luminosity profiles. An- 
dredakis et al. (1995) found a correlation between the value 
of n and the morphological type of the galaxies, in the sense 
that early types show larger values of n than do late-type 
galaxies. Exponential profiles have been used extensively in 
order to fit the surface brightness profiles of discs of spiral 
galaxies. Sersic profiles have been also used in other types 
of galaxies, for example, Davies et al. (1988) propose that 
dwarf ellipticals are well fitted with a Sersic law with n < 2. 
Ground-based astronomical images are always affected 
by atmospheric blurring. Many papers have been written on 
the subject and much effort has been invested in the con- 
struction of new optics to minimize its effects. Seeing scat- 
ters light from the inner, centrally concentrated core to the 
outer, more diffuse regions of galaxies, producing a mean 
surface brightness lower than the true values and larger ef- 
fective radii. These effects can change the results of the pho- 
tometric parameters obtained from the fits of the observed 
surface brightness profiles; moreover, the dynamical prop- 
erties inferred from these parameters will also be wrong. 
Although seeing affects all points in a galaxy, its effects are 
more important in the central regions. Seeing effects were ex- 
tensively studied in the case of elliptical galaxies with r^' ■* 
profiles (Franx, lUingwort & Heckman 1989; Saglia et al. 
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1993). These authors showed that the effects of seeing on 
the photometric properties of eUiptical galaxies can extend 
as far as 5 seeing discs. 

Saglia et al. (1993) also showed that seeing effects be- 
come important for distant ellipticals [cz > 8000 km s~^). 
They found that uncorrected fundamental plane distances 
are systematically too high if seeing effects are not taken 
into account. This is due to the small angular size of the 
objects at high distances: for a flat universe with Hq = 75 
km s^^ Mpc~^, an object of 30 kpc has an angular size of 
15" at z = 0.1 and 3" a.t z — 0.5. This means that a typical 
seeing of 1" is equivalent to 1/15 of the size of the object 
a.t z = 0.1 but is 1/3 of the object a.t z — 0.5. The seeing 
will therefore produce important effects on these objects at 
large radii. Thus, the study of the effects of seeing on sur- 
face brightness proflles must be taken into account when the 
photometric parameters of these galaxies are obtained from 
the decomposition of their surface brightness profiles. The 
new generation of ground-based telescopes and the study of 
galaxies at higher redshifts make these kinds of studies very 
important. 

In this paper we present an analytical treatment of see- 
ing effects on Sersic profiles taking into account the elliptic- 
ity of the objects. The mathematical treatment of seeing will 
be given in Section 2. In Section 3 we present the effect of 
the seeing on the photometric parameters. An easy prescrip- 
tion for correcting the parameters measured from the raw 
profiles is given in Section 4. The analysis of these results 
are given in Section 5. 



Elliptical coordinates (5, 8) are the most appropriate for 
our problem and are defined as: 



y 



^cosS 

5(1 -e) sin ( 



(2) 



In this coordinate system, the surface brightness dis- 
tribution, /(r), of an elliptical source depends only on 5: 
/(r) — 1(5). The convolution equation that represents the 
effect of seeing on the surface brightness distribution is given 
by: 

/c(5,e) = (l-e) / C'de' / de'PSF(5',e',e,e)/(c'),(3) 
Jo Jo 

where PSF(5 ,6 ,^,6) is the Gaussian PSF given by: 



psF{^,e ,^,e) = 



27ra2 



(4) 



2.1 Analytical convolution of Sersic profiles 

In the particular case of Sersic profiles, the surface brightness 
distribution is given by: 



7(C) = /(0)e~(4)", 



(5) 



2 MATHEMATICAL ANALYSIS 

The blurring of images by the atmosphere and imperfec- 
tions in telescope optics (seeing) degrades measurements of 
the surface brightnesses of galaxies. The seeing is character- 
ized by the point spread function (PSF) . The PSF gives the 
probability that a photon will hit the imaging device at a 
point different from where it would have hit in the absence 
of seeing. This can be determined observationally by study- 
ing the scattering of stellar light. PSFs are well described 
by Gaussian functions, Gaussian functions with exponential 
wings, linear superpositions of Gaussian functions, Moffat 
functions, etc. (e.g. Moffat 1969; King 1971; Schweizer 1979). 
Among these analytical approximations of PSFs, the most 
widely used is the single Gaussian. The main goal of the 
present paper is evaluate the effects of seeing on the param- 
eters of Sersic profiles in a completely analytical form and 
to take into account the ellipticity of the surfaces brightness 
distribution in this treatment. We develop our analysis us- 
ing a Gaussian PSF and in Section 5 we shall compare this 
with a different PSF. 

Assume a circular Gaussian function of dispersion a to 
model the point spread function: 

Consider a case where, in the absence of seeing, the surfaces 
brightness distribution J(r) of the galaxy is elliptically sym- 
metric. This means that the isophotes of the object all have 
the same constant ellipticity e (e = 1 — b/a, where a and 
b are the semi-major and semi-minor axes , respectively, of 
the isophote). 



where 1(0) is the central intensity and ro is the scale length 
of the profile. Over the major axis of the object, 9 = 0, the 
analytical solution of Eq. (3) for this type of profile can be 
written as: 



«^.»)^>-.)-'*-i:^f^)'i:^<^ 



1,. 2,1 

e — e 



^^^{^ + ^ + t)M{l + l + ^,l + l,Ui 



1=0 

2^ 



(6) 



where M{jj,, v, z) are the confluent hypergeometric functions 
(Abramowitz & Stegun 1964, p. 504). This expression sim- 
plifles if the object is circular: 



m) = /(o)e- 
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(7) 



In the limiting cases, ,f ^ cx) or ct — > 0, the asymp- 
totic expression of the confluent hypergeometric function 
(Abramowitz & Stegun 1964, p. 504) can be used to recover 
the original Sersic expression: 



/.(e,0)=/(0)e <4)' 



l + O 




(8) 



where (^^/2(7^)~^ quantifies the differences between the 
Sersic profiles unaffected and affected by seeing in those 
asymptotic limits. 
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Figure 1. Effects of seeing on the ellipticity. A model with Te/cr ■■ 
2 is shown. 



Figure 2. Effects of the seeing on the central intensity /c(0) for 
different values of n. Three different ellipticities are shown, e = 
(full line), e = 0.25 (dashed line) and e = 0.5 (dotted dashed line). 



2.2 The effect of seeing on the ellipticity of the 
isophotes 

In the absence of seeing, by construction, all isophotes of the 
Sersic profile have the same ellipticity, whereas the presence 
of seeing tends to make them circular. Using the isophote 
condition, I{£,,0) ~ I{£,,tv/2) — the expression over the mi- 
nor axis is written on appendix A — it is possible to derive 
an implicit equation that gives the variation of the ellipticity 
with the radial distance: 



<0 



-2( |) ln/(e,a^n,ro,C,e(C)) 



(9) 



where f is given by: 

/(e,o-^7i,ro,^,e(0)) 



2 '.(T ' 



2x;r(;+i) 



Z^fe=o fc! ^ ro ' Z^;=oi!'^ ' r(i+i) 

lA^n = nl. (l-e)" H(l-e(C))' "^ V" ^ 2 ^ 



roc 

Jo 



-(|^)".'i+i„-i(^^^)V„(iii^iii)>Lii^)]~i (10) 



'C 



where I„{x) are the modified Bessel functions (Abramowitz 
& Stegun 1964, p. 376). 

Figure 1 shows the radial variation of the ellipticity due 
to the seeing. Note how the seeing affects the central points, 
rounding the isophotes, whereas its effects are progressively 
less in the outer regions of the profile. 

The equations that we have presented for the Sersic 
profiles can be immediately generalized to almost all the 
experimental profiles (see the theorem in Appendix B). 



3.1 The effect of seeing on the central intensity 

To study this effect we use Eq. (6) and apply the fact that 
at ^ = the confluent hypergeometric function satisfies 
M{^, v; 0) = 1. This expression can then be written as: 



/c(0) = /(0)(1 



fc=0 



k\ 



f V2a 



ro 



1 + A 

2n 



a - 

^2' - 



— ■ 1- 2f 



(11) 



where 2F\{a,h;c,z) is the hypergeometric function 
(Abramowitz & Stegun 1964, p. 556). 

As shown in Figure 2, the central intensity of the profile 
affected by seeing decreases monotonically when the seeing 
size increases. The central intensity of profiles with larger 
values of n decreases more rapidly than for low n as is ex- 
pected because of the higher central concentration of these 
profiles. 

The central concentration of the object is also depen- 
dent on the ellipticity. In Figure 2 this relation is also shown. 
Larger ellipticities are more affected by seeing. 



3.2 The effect of seeing on the effective radius 

The Sersic profile can be written in terms of the effective 
radius and the effective intensity: 



/(r) 



/elO 



{^y 



(12) 



The constant bn is chosen such that half the total lu- 
minosity predicted by the law comes from r < Ve- &n can be 
well approximated by the relation 6„ = 0.868n — 0.142. /e is 
the intensity at the effective radius. 

The relation between Je, re and /(O), ro is given by: 



3 THE EFFECTS OF SEEING ON THE SERSIC 
PROFILE PARAMETERS 



m 

and 



/elO" 



(13) 
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Figure 3. Effects of seeing on the effective radius, r^. Three 
different ellipticities for the source are shown: e = (full line), 
e = 0.25 (dotted line) and e = 0.5 (dashed line). 



Figure 4. Values of the parameter rj at rj as function of the 
ratio r^/cr for different values of n and ellipticities: e = (full 
line), e = 0.25 (dashed line) and e = 0.5 (dotted-dashed line). 



ro 



(fen In 10)" 



(14) 



The seeing effect on effective radius can be obtained 
from the conservation of luminosity by the convolution; 



L^irl) = i(r,), 



(15) 



where L{re) is the luminosity of the source inside re and 
L'^{r'e) is the luminosity obtained from the object affected 
by seeing, measured inside its effective radius. For a circular 
object, e = 0, we have 

L^irl) = 2ti j rh{r)dr (16) 

and L(re) = I{Q)iTrl[nr{2n)]/[{br, InlO)^"] for a Sersic pro- 
file. Equation (15) can then be written analytically for a 
circular system as the implicit equation; 



nT{2n) 



:^HV^- 



/.—() 



k\ 



^(fe„lnlO)2" 

^^ (-)' 1 (-^), (rf\ 

Z-^i^o n i+i (i)i v^CT^y 



ro 



r 1 



2n 



(17) 



where (a); is the Pochhammer symbol: (a); = r{a+l)/T{a). 
Figure 3 shows that the effect of seeing is to increase the 
effective radius. This effect becomes more important as n 
increases. The ellipticity effect is also shown; however, for 
e 7^ there is no easy an analytical form, so the results with 
t ^ that are shown in Figure 3 were obtained numerically. 
Greater ellipticities imply greater effective radii, these dif- 
ferences are more important for greater values of n. This 
result is as expected due to the diminution of the central 
intensity by the ellipticity effect. 

It must be noted that our measurement of effective ra- 
dius has been obtained over the semi-major axis. Some au- 
thors use as radial distance the magnitude r* — vafe, in this 
case, the effective r adius of th e object affected by seeing is 
given by r" — r^-v/l — e(re), where e(re) can be obtained 
from Equation (9). 



3.3 The effect of seeing on the parameter n 

To quantify the effect of seeing on the parameter n we use 
the parameter; 



Vi^) 



1/(0 1^/(0 



di 



1(0)- 



(18) 



This parameter is defined in such a way that i]{(^) — n for 
all values of ^ if /(^) is a Sersic profile (Eq. 5). So ri{^) is 
equivalent, locally, to the parameter n of the Sersic profile. 
Figure 4 summarizes the values of this parameter at r^ for 
different n values. It is easy to see that this parameter is 
the most affected by the seeing. Indeed, 77(0) — 0.5 for any 
profile affected by a Gaussian seeing. It should be noted that 
Sersic profiles with n = 0.5 are Gaussian profiles and then 
its convolution with a Gaussian PSF gives another Gaus- 
sian. The principal conclusion is that seeing effects always 
produce a surface brightness profile with a smaller value of n 
than the actual value. Again, as n increases, the parameter 
is more affected. The effect of the ellipticity is to decrease 
the value of n but the changes are not so important. 



4 A PRESCRIPTION FOR SEEING 
CORRECTIONS 

Here we present an easy prescription based on the use of 
the plots of Figures 2, 3, 4 and 5 (see below). This proce- 
dure permits the parameters of the Sersic profile (seeing-free 
quantities) to be obtained using the observational surface 
brightness profile. In summary, observers should; 

a) Determine the FWHM (Full Width at Half Maximum) 
of stars by fitting a Gaussian, a is related to the FWHM by 
a = FWHM/V8]n2- 

b) Measure r^ along the semi-major axis solving the im- 
plicit equation L'^(r'i) = (l/2)L'^(oo). This can be done 
(without any assumptions) directly from the raw images. 

c) Determine r]{rl) numerically using the expression; 
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Figure 5. Values of the ellipticity of the isophotes at r^ as func- 
tion of the ratio r^/cr for different values of n. 



Vir-e) = 



1 Jc(r-g) 



In 



/c(0) 



(19) 



d) Evaluate the value of n which corresponds to the point 
(»7(re), rl/a) using Figure 4. Suppose, as a first approxima- 
tion, that the value of e corresponds to the value of e{r'^)- 
Note that r;(^) is the parameter less affected by the value of 
e, so the approximation is good. 

e) Recalculate the value of e more accurately using the 
value of n obtained and Figure 5. Figure 5 shows the values 
of e{r1) for different values of n. 

f) Obtain the value of r^ using Figure 3. 

g) Obtain the value of 1(0) using Figure 2. 

Observers wishing to be more precise can use the formulae 
instead the figures. 



5 DISCUSSION 

The study presented here has assumed a Gaussian PSF for 
the seeing. The real observed PSF is not exactly Gaussian. 
The theory of atmospheric turbulence predicts the PSF to be 
the Fourier transform of e'^'^*''^'''^ (Fried 1966; Woolf 1982), 
where 6 is a scaling parameter. Saglia et al. (1993) general- 
ized this result, they assume a PSF that is the Fourier trans- 
form of e~ ''''•' . The Gaussian PSF is a particular case with 
7 = 2. Their observational PSFs were in agreement with 
the theoretical PSF inferred by the turbulence theory. They 
obtain a Gaussian FWHM 4.67% greater than the FWHM 
of the turbulence PSF with 7 = 5/3. This systematic error 
will be transmitted into the parameters le, r^ and n. We 
have computed these parameters, varying a in our PSF by 
4.67%, and have found that Je is underestimated by 7%, Ve 
is overestimated by 7% and n is underestimated by 4% in a 
systematic manner with respect to the initial values. Due to 
the systematic character of these errors they can be easily 
taken into account. 

The parameter most affected by the seeing is n. This 
has important consequences because n serves to classify the 



Figure 6. Normalized luminosity den- 

sity J{(^) = JiQrJ" /[f^^^I{0)] vs. the radial coordinate C/ve 
for Sersic profiles with n < I: n = 0.2, 0.3 and 0.4 (dashed lines), 
n = 0.5 (solid line) and n = 0.7,0.8 and 1.0 (dotted lines). 



type of profile and is related to the central luminosity con- 
centration. Graham et al. (1996) and Jerjen, Bingelli & Free- 
man (2000) found a correlation between n and galactic type. 
Thus, dwarfs ellipticals show the smallest values of n and 
cD galaxies have the largest values. Between these extremes 
are located the ellipticals and the bulges of spirals. The pa- 
rameter ri{^) defined in Section 3.2 gives information locally 
about the value of n over the profile. Figure 4 shows rj{r^) 
for different values of the seeing. It can be observed that the 
seeing always produces rj{r^) < n for all the values of the 
r'^/a ratio. If seeing is not taken into account, the value of 
n that can be measured from the profiles is always smaller 
than the real one. Usually, fitting procedures avoid the cen- 
tral points in order to remove seeing effects on the profiles. 
This is clearly not sufficient to recover the real value of n 
(see Figure 4). 

One physical restriction to the values of n is given by 
the luminosity density. For a homologous triaxial ellipsoid, 
the luminosity density (Stark 1977) associated with a Sersic 
profile is given by: 



i(C) 



f2 KI{0) 



1/n 



(*)%* 



(■-'(,('-(') 'Jf,(20) 



where K = bnln 10 and i^'^ is a constant that depends on 
the 3D spatial orientation of the object. We have calculated 
the luminosity density for Sersic profiles with n < 1 (see 
Figure 6). For n < 0.5 the density has a depression in its 
central parts. This represents an unlikely physical situation. 
Nevertheless, the seeing effects prevent the measurement of 
n < 0.5 for objects with n > 0.5 (see Figure 4). 

The effects of seeing on the parameters le , r^ and n de- 
pend on the intrinsic ellipticity (e) of the source. Thus, as- 
suming that the object has e = 0, when it really is elliptical- 
symmetric, results in the central intensity and n begin un- 
derestimated whereas re is overestimated. These effects are 
not negligible and they are more important when e increases. 
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6 CONCLUSIONS 

We have developed an analytical study of the seeing effects 
on Sersic profiles. The seeing PSF was modelled by a Gaus- 
sian function, fn this analysis we have taken into account 
the intrinsic ellipticity of the objects. Our main results are: 

(i) The convolved surface brightness profile along the ma- 
jor axis of the object can be expressed as a double series of 
confluent hypergeometric functions. This result is very gen- 
eral and can be applied to nearly all the experimental surface 
brightness profiles. 

(ii) We have obtained an implicit equation to evaluate 
the effect of seeing on the ellipticity of the isophotes. The 
rounding of the isophotes depends on n,re,(y and e in a 
unique way. 

(iii) The parameter most affected by the seeing effect is n. 
The observed Sersic profiles show smaller values of n, due to 
the seeing effect, than the real ones. Greater values of n are 
the most affected. Also, for n < 0.5, the luminosity density 
associated with a Sersic profile has a depression in its central 
parts. This represents an unlikely physical situation. 

(iv) The seeing effects on the parameters of the Sersic pro- 
file depend on the intrinsic ellipticity of the object, therefore 
it is necessary to include it when seeing effects are studied. 

The results described here clearly show that seeing effects 
are important when one tries to measure accurate values 
of the parameters of a profile affected by the seeing. These 
results have to be taken into account for sources with a 
low re/a relation as expected for medium- and high-redshift 
objects. 
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APPENDIX A: ANALYTICAL EXPRESSION 
ON MINOR AXIS 

Using Eq. (3) we can also get the equation for the minor 
axis, 9 = 7r/2, 






\/2cr 



ro(l-e) 



1=0 n 



WW 



' £(f±i2r a + 1 + il 

V{1+1) *- y^ ^^ 2n) 



M(/ + l+A,Z + l,i[^^]^ 



(Al) 



However, this expression is 

divergent for |(e^ - 2e)/(l - ef\ > 1, so that for e > 0.3 
the use of the integral expression is required: 



where F2(^,5 ,<^,^) is given by 



1 (i-o^s 
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(A2) 
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In the limiting cases, ^ — > oo or a ^ 0, for e < 0.3, 
the asymptotic expression of the confluent hypergeometric 
function (Abramowitz & Stegun 1964) can be used again to 
recover the original Sersic expression: 



(f,f)^/(0). 



l + O 






.(A4) 



APPENDIX B: THEOREM ON GAUSSIAN 
SEEING 

Assume a surface brightness distribution with elliptical sym- 
metry, f (5) . If this distribution can be written as the power 
series 



/(c) = /(o)^«(fc)e 



/3fc 



(Bl) 



with a region of convergence equal to < .f < oo , then the 
analytical solution of the convolution of /(^) with a Gaussian 
PSF over the major axis is: 



© 1994 RAS, MNRAS 000,|i|-| 



The effects of seeing on Sersic profiles 7 



^Md 



kfiy M-^ ,,.,,- -.Ilk 



k = 



2\' r(;+i), 



3k^ 
2 j 



M(/ + l + f,/ + l,i(|)'). (B2) 

For e=0, the convolution can be written as 



*,— 



M(l + f,l,i(i)^). (B3) 

For the Sersic profile, a{k) = [(-)''7A:!](l/ro)''''", /3 = 1/n. 

Examples of profiles where the theorem is not applica- 
ble are the Hubble profile (Hubble 1930) and the Freeman 
bar profile (Freeman 1966). The first one does not admit a 
convergent power series over the entire interval of definition, 
and the second one has a point of non-differentiability that 
avoids a power series expansion. 

This paper has been produced using the Royal Astronomical 
Society /Blackwell Science I^T^ style file. 
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